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A method is presented for obtaining the backscattered radiation intensity from an 
idealized target illuminated by a beamwidth limited source using the concept of differential 
reflectivity. In general, this concept provides a method for determining t lie fields reflected 
from a body illuminated by an arbitrary source or antenna pattern when the reflection co- 
efficients for a plane wave incident on the body are known or approximated. Fresnel's 
reflection coefficients are claimed to be appropriate approximations for the class of convex 
targets characterized by either a large radius of curvature or consisting of lossy material. 
In the application of this method to the case of a, large spherical body, considerable computa- 
tional simplification results. Application of results to appropriate experimental data yields 
a minimum value of 1.53 for the average dielectric constant of the moon's surface material. 
Extension of the method to statistical problems is indicated. 



1. Introduction 

In recent years numerous attempts have been made to obtain information on the proper- 
ties of such distant bodies as the moon and neighboring planets [Beckmann, 1963]. A major 
method involves analyzing and interpreting the backscattered radar return from these bodies. 
This backscattered radar return may be analyzed from a number of different aspects, depending 
on the nature of information sought. Through statistical analysis, such information as pulse 
ensemble averages for discrete times during the return as well as autocorrelations, cross correla- 
tions, and variances thereof may be found. But even for a body as close as the moon, the 
information that can be obtained from such analyses will represent gross averages over large 
areas of the body. 

Information of particular interest is the value of each of the electromagnetic properties, 
such as e, /jl, and a, which even in simple nonstatistical problems are extremely difficult to 
separate from each other. However, this separation can be achieved when the data is examined 
in the light of some simplifying theory capable of predicting the return for sharply defined 
ideal conditions. An analysis of the return data together with the theoretical solution should 
then yield average values for the electromagnetic quantities sought. 

This paper presents a solution describing the radar return from a simple geometric shape 
and an application of this solution to the case of an idealized moon. In the case of the moon, 
the idealization is necessary to de-emphasize the contribution of factors such as roughness and 
local variations in electromagnetic properties. The derivation of the integral representing the 
radar return is rigorotis, but certain simplifying assumptions are necessary to carry out the 
integration in closed form. 



1 This work sponsored by the National Aeronautics and Space Administration under Grant NsG 129-61. 
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2. Theory 

2.1. Concept of Differential Reflectivity 

A concept which will prove useful in the development of the theory outlined in this paper 
is that of "differential reflect ivity." This is a dyadic quantity. When this dyadic is multi- 
plied by a differential surface area and the vector field incident on that area, there results an 
expression of the contribution of that surface element to the scattered field at an arbitrary 
observation point. The differential reflectivity is a function of the following: 

(1) Location and orientation of the surface with respect to the observation point. 

(2) Properties of the two media separated by the surface. 

(3) Frequency of the incident radiation, co. 

For the case of steady state incident radiation the reflected Hertzian potential field, for instance, 
may be described by 



U r (r l} co, 0= I I tfOi, r , <a)'TLt(r Qf co, t)dS 

So 



(1) 



where r x and r denote the coordinates of the field point and the surface points, respectively; 

-» -> -> -> 

n ( (r , co, t) is the incident vector field at the surface S ; and a(r u r , w) is the dyadic differential 

reflectivity. 

Consider now that the radiation incident on a surface is from a pulsed, "conical" source. 
A "conical" source is understood to be one whose radiation is limited to a cone of vertex angle 
a. Within the cone the electric and magnetic fields are uniform over any spherical surface 
centered on the source. Exterior to the cone they are identically zero. Let the source radiate 
between times /=0 and t=T a sinusoidally varying carrier wave of angular frequency co . 
At some later time the outwardly traveling wave impinges on an infinite plane surface (z=0 
in fig. 1) which separates all of space into two semi-infinite regions, each of which is filled with 
a homogeneous medium. The portion of the plane intercepted by the cone is given by S . 
The intensity of the backscattered field at an arbitrary observation point is desired. 




X direction 



Figure 1. Reflection geometry. 
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If one replaces the incident wave packet with an infinite set of steady state (conically 
bounded) incident waves obtained by means of the Fourier transform, one may solve for the 
contribution of a typical member of the family of steady state waves and, finally, sum up or 
integrate these contributions over the whole family to obtain the solution to the original 
problem. It may be observed that WeyPs method of expansion into plane waves [Stratton, 
1941, pp. 577-582] is not applicable to the beamwidth limited case. 

Using the expression for the differential reflectivity, the contribution to the total field at the 
observation point due to a conical bounded typical member of the family of steady state waves 
can be written as 

-» -> C C K "^ ~* "* "* 

ML r (r u co , «)= ff(r u r , co) -n t (r , co , w)dS (2) 

where o-^, r , co) is the differential reflectivity and II i(r , co , co) is the complex phasor for the inci- 
dent field due to the steady state wave of angular frequency w. The time dependent total 
scattered field is then 

n r (r lt u 0l =— j °° e~ i<at J a(r h r 0l co) -II, (r , co , co)dS J A*. (3) 

So 

2.2. Theory of Differential Reflectivity 

With attention focused on only one member of the family of steady state waves, a deriva- 
tion of the theory involving the concept of differential reflectivity will now be shown. Let the 
Hertz vector due to the component steady state incident wave be described by 

II,(/ C . u,t)=a,C ^-e-" t (4) 

— f 

for all points r on the surface S . On the remainder of the infinite plane surface it is identically 

zero. Here a T is the unit vector in the II t direction, Co=Co(«) and relates to the source 

strength, k 2 =k 2 (oS) is the propagation constant in the incident medium, and R =r —r s , where 

-> 

r s is the radius vector from the origin to the source point. 

If the origin is taken in the infinite plane surface (the x-y plane) of which S is a region, and 
the source has the rectangular coordinates (0, 0, z s ) as shown in figure 1, so that 



there can be written 



n t (r , co, O = 0»£V-' W< I -g7- &(xo—Xo)Hyo—yo)dxody' (5) 



So 

where 



Using the Fourier integral expansion one may write 

CO 

d(xo-x' )5(y -y'o)=^ \ \ e tu < x *~ l b +"**<' duds (6) 

— 00 

so that the incident field at each point, r , of the surface is given by 

CO 

n,(r , a,. 0=^r f T 1 !? [ 0* ^ c o ei " (xo " J » )+ "' < " )_; '" > ^'^] dx '^h- (7) 
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If now r=^ i xl-\-yl J rzl there can be obtained by analytic continuation the expression 

CO 

n,(r, », 0=^2? f f TgT [ fp^o^^^^l cfetfi/i (8) 



So 

where 



Wi=u(x — x' o )+v(y — ?/o) — ^oV^2— (^ 2 +^ 2 ). 
-> -» 
The expression (8) which reduces to (7) for r=r is subject to the following interpretation: it 
can be considered to be the field due to an infinite collection of plane waves, symmetrically 
distributed about the local normal to the incident wave front, which combine at a point on the 
reflecting surface to yield the net incident field due to the original source. The propagation 
constant associated with each of the plane waves is so determined that the entire collection adds 
to a two-dimensional delta function at the point in question. 

Consider now a Hertzian plane wave having a propagation vector with components 

[u, v, ^kl—(u 2j rv 2 )] and polarization in the direction a r to be reflected from the surface. The 

reflected plane wave will be described through the use of a dyadic reflection coefficient 

a s <- 

V (u, v) . At an observation point Q( r i) the total reflected field due to illumination of the surface 

S by the infinite set of plane waves will be given by 

n r (r„ co, t)= |J ([^ jjv(u,v)e iW rdudv^^C ^p e-^dx'dVo (9) 

So —<» 

where 



Comparison with (1) yields the differential reflectivity 



00 

— 2 ff V(u,v)e iW rdudv. (10) 



° 4 

The form of the components of the reflection coefficient V will depend on the nature of the 
surface S . If the surface is spherical they may be derived from Mie's solution for a plane wave 
incident on a sphere [Stratton, 1941, pp. 563-567]. For an infinite plane surface they reduce 
to the ordinary Fresnel reflection coefficients. 

2.3. Derivation of the Components of a 
For purposes of computation it will prove useful to evaluate the vector quantity 

2=^2 ff V(u,v)-aX\e iW rdudv=a.11 7r C . (11) 

— CO 

Referring to figure 1, let there be defined two coordinate systems having their origins at the 
point P on the reflecting surface. The Q system will be defined by the orthogonal set of unit 

vectors a t , a r and a n and the K system by the orthogonal set a P , a N and a n . Here a n is the 

-> -> 

positive unit normal to the reflecting surface, a r is in the direction of the projection of B Y on 

->-»-> -> 

the tangent plane through P, and a t =a r Xa n . The vector a P is in the direction opposite to 

-> ->->-> 

that of the projection of k 2 on the tangent plane and a N =a n Xdp- 

These two reference systems are in addition to the primary reference system in which, in 
a generalized problem, the surface normal changes direction as one traverses the surface under 
consideration. 
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The vector (11) may be resolved as follows: 

00 00 

X} = t~ 1 a P VppC P e iW r dudv+ -r— 2 a N V NN C N e iW r dudv 

+ ~h ((a n V nn C n e iW rdvdv=£+T>+lt (12) 

where (M7o=flpCp+G;vC2v+ a »C^- For plane Il-waves reflecting from an infinite plane surface 
it can be shown that 



V pp (a) = — V nn (a) = — 



2 COS a— -yjn 2 — sill 2 a 



( — m 2 cos a 



+V^ sm 



( — ) cos a— V^ 



A 



7 \ (13) 

( — ) cos a-\--\'n 2 —s'm 2 a 

\M2/ 

where a=angle of incidence = sin"" 1 ! — p — )> and n=index of refraction=T- In this case V, 

and F iViV are identical with the Fresnel reflection coefficients. Off-diagonal components of V 
are identically zero by virtue of the choice of reference system. 
Let 

u=\ cos P=k 2 sin a cos (3 

v=\ sin t3=k 2 sin a sin 

?/i — 2/o = Pi sin ft 

Pi=Ri sin #! 

g^BiCOS^. (14) 

On substitution into (12) one obtains 

s=«« ^2 r F,«[a(x)] r r * ^ c ° s ^i-»+'*i*« °™ «^i x^x 

=o„ ^ J F M [a(X)]-/o(Xpi)e'^ cos "XtfX. (15) 

Also one can obtain 

=a t £ J*" [/ l (X)«7 (XPi) +./ 2 (X) J 2 (X Pl ) ]e'^ cos "XrfX 

+«r^f L/'i(X)J (Xp 1 )-/ 2 (X)J 2 (X Pl )]«^ cos »XrfX (16) 

where 

./ 1 (X)=Fppla(X)]+F^ v [«(X)J 
/ 2 (X) = F PP [a(X)]-F^[a(X)| 
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C t = (a t • a r ) C = (a t • a P ) C P + (a t • a N ) C N 

C r = (a r • a*) C = (a r • a P ) C P + (a r • a N ) C N . 

-> -» 

x\pproximate evaluation of 2 n and 2 r as given in (15) and (16) may be made by use of 

the saddle-point method [Brekhovskikh ; 1960] which, under the assumption k 2 Ri sin 2 0!>>1, 

yields 

^ T =at £ y ™ ( * i} ^r cos ^ +a ' £ ypp( ^ i} ^r eos ^* (17) 



From the expression S^c-a^Co we now find 

ko £ ili 2 R \ 

<Ttt= 2TTl VnN ^ ~#~ C ° S $1 

h ~ T7 m \ eik2Rl a 

I, pi^2 R l 

anr>= ~Si Vpp{dl) ~eT cos ° l (18) 

as the components of the differential reflectivity with respect to the designated coordinate 
system. 

Similar results may be obtained for surfaces of arbitrary curvature by use of the appro- 
priate form of the reflection coefficient. It is realized that evaluation of the reflection coefficient 

V for a plane wave incident on an arbitrary surface may be a problem in itself. 

2.4. Discussion on a 

The differential cross section derived above is readily applicable to the problem of scat- 

tering from an infinite plane surface. It is noted from the form of the vector 2 that the field 
strength at right angles to the surface normal is zero due to the cos X factor. In line with 
our original assumptions then, a plane surface can be considered as a collection (infinite) of 
independent secondary sources without mutual coupling. This is a remarkable result inas- 
much as seldom in 3-dimensional electromagnetic scattering problems can neighbor-to-neighbor 
interactions be ignored. If, now, a portion of the surface illuminated has the cross-sectional 
profile shown in figure 2, it is apparent that there will be a contribution to the fields at both 
B and C due to the secondary source at A and vice versa. The strength of the secondary 
sources everywhere on the surface will thus be affected. The possibility of secondary reflec- 
tions such as occur at B may be negated by requiring that the bodies illuminated be convex. 
The situation at C cannot be so simply disposed of. However, it can be seen that if the 
radius of curvature of the surface is large, the contribution at C due to the source at A will 
be small. In addition, if the material composing the target body is lossy, the attenuation 
suffered in traversing the path A to C will be so high that the re-radiation from the C surface 
element due to the illumination of the A surface element will be negligible. The method 
outlined herein is, therefore, an excellent approximation for the solution of scattering problems 
involving large convex bodies or any highly lossy convex body. 
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Figure 2. Reflection from nonplanar surface. 




Figure 3. Geometry for reflection from a sphere. 



3. Application to the Idealized Moon 

An application of the method to the case of a large spherical object such as the moon can 
easily be made with the aid of a few approximating assumptions. In particular, the spherical 
object will be assumed to be a smooth sphere of large radius, a, and composed of a homogeneous 
lossless substance. Geometry pertinent to the problem is shown in figure 3. The origin of the 
spherical coordinate system is located at the center of the sphere with the source and the re- 
ceiver being located at r=D and 0=0. Additionally, let 6=^ and <p=0 define the positive 

x-axis; and 6=0 define the positive 2-axis. 

Let the source be an electric dipole oriented in the ^-direction so that 



Ui=a x C 



e ik 2 R 

~H7 



(19) 



where B Q is the distance from the source to a point on the surface of the sphere. It is readily 
seen that 



C t = — C sin <p 
C r = — C cos <p cos 6 
C n =C cos <p sin 6 



(20) 



Realizing that the proper reflection coefficients to use in this case would be those obtained 
from Mie's solution to the problem of a plane wave incident on a sphere, the assumption is now 
made that due to the very large radius of the sphere compared to the wavelength of the incident 
radiation, these reflection coefficients may be replaced by their limiting values as the radius of 
the sphere increases without limit. These limiting values are, ipso facto, the Fresnel reflection 
coefficients as given in (13). 

From (1), (18), and (20) there is obtained 

a 2 fa C T ~* "* 

II r = — (7 ^— - cos O^dt sin (pV NN (Si) J ra T cos <p cos 6V PP (6i) 

8 

+a re cos (p sin dV PP (6i)] — ^ ^— sin BdBdcp (21) 

llQ ill 
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where E x is as before, the distance from a point on the surface to the observation point. If in 
this expression one replaces B x by R so as to obtain the backscattered field at the position of 

the source one finds, on performing the integration over <p, that (21) vanishes. However, as 

-> -» 

will be shown, the E and H fields at this point do not vanish. We will then be able to compute 
the net power incident on an antenna of given effective area. 

-> -> 

The elementary far field contributions to the E and H fields at the observation point are 
given by 

82?(Bi)= -lda B X (a R X8U r ) 

bH{R x ) = uk 2 e 2 (a R X 8IL r ) (22) 



where 



and 



a R =a r sin 6}-\-a n cos 6 X 



8U r =—C Q 7z — r cos B x [a t sin ^PWfflO+a,- cos p cos 0Fpp(#i) 

+a n cos <? sin ^ PP (^)] ^— -^— sin BdSdp. (23) 

One obtains, on making the indicated substitutions and integrating over <p, the remaining 
components of E and 7? 

E x {R,)=^p\ D2 sin B x cos ^[FppWO cos 2B l +V NN {B l )]dB l (24) 

ff/fl )= , 2 / 2<0 n2 sin 9, cos ^[Fpp(^) cos 2^+^(001^1 (25) 



and 



where R ~D and 0i ~0 have been used as approximations since Z»>a. 

In order to obtain numerical results for the problem at hand, the Fresnel reflection coeffi- 
cients V PP and V NN are replaced at this point by the expressions 

Vp P =L 3 cos 2 6>+Z 4 cos 0+1 

V'm^Lx cos 2 0+Z 2 cos 6— 1 (26) 

where 

Z 1 =~0.7 

£ 2 =P+1.7 
i — *-+± 



L,—[L 3 h+l] 



l+n 
and 

A = cos & 
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P=F™(0)=i-r^if Mi=/i 2 (27) 



Figure 4. Comparison of exact and approximate 
reflection coefficients. 




where 6 b is that angle satisfying the relation V PP (6 b ) = 0. The approximation is based on match- 
ing the curves V PP (6) versus 6 and V NN (d) versus 6 at three points, namely, at 0=0°, 0=90°, 
and d=6 b . Figure 4 shows the approximate reflection coefficients V PP and V NN in relation to 
the exact coefficients for the case: 



Let now 



jUi = jLt 2 , o-i = a 2 = and n =1.5. 



a 2 n n b s 

**~ 2D* 



II, 



2J9 2 



(28) 



where 



r=-d 

Jo 



<? m 2 (D-«cosa 1 ) sin ^ cos ^[^(flj) cos 2<9 1 +T 7 : N ,v(0 1 )k/6> 1 



On substitution of the expressions for V PP and V NN there results, after integration and dropping 
of terms involving powers of lyjTT. higher than the first, 



e l2k 2 L 



2k 2 a 



[ e -i2k 2 a cos e a j(e a ) —f oe -^2a] 



(29) 



where /o=2p and 

f(d a ) = 2L 3 cos 5 e a +2L 4 cos 4 d a +(L l -L s +2) cos 3 d a +(L 2 -L 4 ) cos 2 0,-2 cos 6 a . (30) 

From the relation 



is obtained 



S=iRe(£X/i*) 
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(31) 
(32) 



where F(d a )=f 2 (6 a ) J rj : l—2f f(d a ) cos [2k 2 a(l— cos 6 a )] and # is the average steady state power 
returned for a fixed beamwidth 26 2a of illumination, corresponding to an illuminated area of 
the spherical surface subtended by the polar angle 6 a (measured at the origin). The angles 
d a and B 2a are shown in figure 5. 

Now it is desired to evaluate C\ in terms of normal radar quantities in order to obtain the 
power received P r by the radar. If P T is the peak power radiated by the antenna, G T is the 
gain of the transmitting antenna over an isotropic antenna, and G SD is the gain of a short dipole 
over an isotropic antenna (—3/2), then the power which must be radiated by the short dipole 
(W S d) to yield the same power density in the main lobe is given by 



Then [Stratton, 1941] 



W s 



CI 



G S d 3 



3W a 



47TCO €2Al2'VA l 262 



(33) 



(34) 



The received power, P r , is given by 



Pr = 



47T 



(35) 



where G R is the receiving antenna gain over an isotropic antenna and X is the wavelength at 
frequency 8. 

Using (32), (33), and (34), (35) becomes 



Pr( 



P T G T G R \ 2 {iva 2 ) F(6 a ) 



(4tt) 3 L> 4 



4 



(36) 



It can be observed that the third term of F(6 a ) as given in (32) is extremely sensitive to 
variations of 6 a . This is because of the large value of k 2 a for a body of the size of the moon 
at microwave frequencies. Figure 6 shows F(6 a ) plotted versus 6 a for three distinct situations: 



cos [2fc 2 a(l — eos0 a )] = l, 0,-1. 



(37) 
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Figure 5. Beamwidth-limited illumination. 
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Figure 6. Relative amplitude of the envelope of 
S versus a . 
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Noting that F(0 a ) converges to only one value for a =9O°, (36) becomes 



^(90°)= P ( ^ff (pVa 8 ). (38) 

The last term in (38) is just the radar scattering cross section of a sphere, that is the power 
reflection coefficient, p 2 , times the geometrical cross section, ira 2 . The first part of (38) is just 
the standard two way radar range equation [Kerr, 1951]. Therefore, the results that are ob- 
tained in this paper using the concept of differential reflectivity are identical with those obtained 
under the same conditions using the Mie solution. However, using the method presented 
in this paper one can obtain the steady state power returned for any specific antenna pattern 
including partial illumination of the moon. This cannot be done by any other known method. 

4. An Application to Experimental Data 

Since it is often possible to obtain information about the transient solution of a problem 
from its steady state solution, an attempt is made here to predict the possible transient solution 
of the problem under consideration without performing the integration indicated in (3). 
Letting 6 a be considered as a function of time it is seen that f in F(6 a ) is independent of time 
and will be accepted by the receiver ;/(0 u ), is a slowly varying function of time and would, no 
doubt, show up in the receiver output. The last term, being a high frequency term (of the 
order of the transmitted frequency), would not be passed through the receiver because of its 
band pass characteristics. The rise time of the first two terms, however, would be too high 
for an ordinary radar receiver, and the initial portion of the return signal would be limited by 
the step response of the receiver and/or the transmitted waveform. 

In addition, in the practical case of return from the moon, roughness over the first few 
Fresnel zones would cause a considerable modification in the received initial slope giving, 
possibly, an ensemble average slope of the received pulses less than that of the receiver step 
response. In any case the maximum of the average returned power in an experiment should 
be indicative of average electromagnetic properties of the moon's surface for the assumption 
that the surface roughness is negligible. If the roughness were not negligible, the maximum 
average value would be less than that for the case involving no roughness. Correspondingly, 
the average electromagnetic properties obtained by matching the maximum ensemble average 
of actual returns with the maximum indicated by this theory (for a smooth moon), would 
yield the minimum average values for the actual moon. 

The received power will now be calculated and applied to the data obtained from a pulse 
radar lunar reflections experiment by Mathis [1963]. By the previous discussion only the 
DC and slowly varying terms of F(9 a ) will be used in the computation of P r . The radar 
parameters of Mathis' experiment are listed in table 1. The maximum value of F(8 a ) is 

U 2 o+f(6 a )] mSiX =8p 2 . (39) 

Then (36) becomes 

Pr m ^°^0 (2pW) (40) 

and on letting 

P T =2.5X10 6 W 
Z?=2.34X10 5 miles 
a=1.08X10 3 miles 
G R =G T =Z7.5 dB 
Pr raax (dBm) = -77.3+20log 10 p. (41) 
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Now using Mathis' [1963] maximum value for -Pr max of —96.8 dBm, one obtains \p\ =0.106 and 
upon letting m = n 2 



l±|p| 

"i-IpI" 



l24 =Vs- 



(42) 



This yields a minimum average value of e for the material composing the moon's surface. 
Thus 6 avg =1.53e , a value which is consistent with those obtained by other investigators. 

Table 1. Trinidad test site radar parameters 
[Mathis, 1963] 



Transmitter 



Receivers 



84 ft parabolic dish on Az-El mount 

37.5 dB gain 

2.25° beamwidth 

Transmit one polarization 

Receive transmitted and orthogonal polarization 

Continental electronics AN/FPT-5 

425.0 Mc/s 

2000 fxsec pulse 

30 or 30.2 pulses/sec 

2.5 MW peak power (nominal) 

Rectangular pulse 

Noise figure 5.5 dB maximum 
Frequency stability— 1 part in 10 in 33 msec 
Gain stability— 0.25 dB in 4 lir 
Dynamic range— 40 dB (linear) 
Bandwidth 4.7 kc/s 



5. Conclusions 

This paper has presented an exact solution to beam-limited reflections of a conically 
bounded spherical wave from a semi-infinite plane through the use of the differential reflectiv- 
ity. In general, the use of differential reflectivity provides a method for decomposing the 
radiation from an arbitrary source into an infinite set of plane waves at each point of the 
illuminated area so that, in order to determine the reflected fields, only the reflection coefficients 
for a plane wave incident on the body need be known or approximated. It has been quali- 
tatively shown that the approximate reflection coefficients for a convex body with either a 
large radius of curvature or consisting of a lossy medium are the Fresnel coefficients. 

The theory has been applied to the case in which a beam-limited spherical wave is incident 
on a smooth sphere. The steady state backscattered power returned from the sphere is obtained. 
The time average steady state power returned displays very rapid fluctuations with the increase 
in the number of Fresnel zones illuminated. Upon matching this expression for average 
power with the appropriate experimental data as obtained from the moon, the minimum value 
for the average dielectric constant is found to be 1.53. 

This method, using the differential reflectivity, proves to be a useful tool which can also 
be applied directly to the case of acoustic scattering. With the use of high-speed digital 
computers, the restrictions imposed herein merely to facilitate computation can be easily 
relaxed. The method may also be easily modified so as to incorporate statistical variations 
in any electrical or geometrical parameters. 
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